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Abstract 



A variety of strong and electroweak interaction properties of the pion and 
the hght scalar a meson are computed in a relativistic quark model. Under 
the assumption that the resulting coupling of these mesons to the constituent 
quarks is identical, the a meson mass is determined as Mg- = 385.4 MeV. We 
discuss in detail the gauging of the non-local meson-quark interaction and 
calculate the electromagnetic form factor of the pion and the form factors of 
the vr^ — > 77 and a — > 77 processes. We obtain explicit expressions for the 
relevant form factors and evaluate the leading and next-to-leading orders for 
large Euclidean photon virtualities. Turning to the decay properties of the a 
we determine the width of the electromagnetic cr ^ 77 transition and discuss 
the strong decay a — > vrvr. In a final step we compute the nonleptonic decays 
D cjvr and B — > an relevant for the possible observation of the a meson. 
All our results are compared to available experimental data and to results of 
other theoretical studies. 

PACS: 12.39.Ki, 13.25.Ft, 13.40.Gp, 13.40.Hq, 14.40.Lb 

Keywords: it and a meson; B and D meson; relativistic quark model; nonleptonic and 
electromagnetic decays. 



I. INTRODUCTION 

During the last few years the physics of isoscalar scalar mesons and in particular of 
the a meson has received an revival of interest due to substantial progress in experimental 
and theoretical activities [1] (for a status report see, for example, Ref. [2]). Originally, 
the light scalar meson a was introduced as the chiral partner of the pion in the two-flavor 
linear cr-model [3,4]. The linear a-model fulfils the chiral commutation relations, contains 
the partial conservation of the axial current (PCAC) and has a manifestly renormalizable 
Lagrangian. In approaches based on the linear realization of chiral symmetry (see, for 
example, [5,6]) the a meson serves to create spontaneous breaking of chiral symmetry, it 
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generates the constituent quark mass and it is a relevant ingredient in explaining low- energy 
phenomenology (tttt scattering, A/ = 1/2 enhancement in K ^ tttt, attractive interaction 
between baryons in nuclear matter, etc.). On the other hand, the use of the linear or 
non-derivative cr-model Lagrangian leads to well-known difficulties. For example, current- 
algebra results cannot be reproduced at the tree level and can only be generated by mixing 
up different orders in the chiral expansion. For this reason, it was suggested [7] that the 
linear o"-model Lagrangian is modified in such a fashion that the results of current algebra 
are already produced at the tree level, while at the same time a clear chiral power counting 
is obtained. This modification is based on a unitary, chiral field-dependent transformation 
of the quark/nucleon field, which eliminates the non-derivative linear coupling of tt and a 
and replaces it by a nonlinear derivative coupling of the chiral rotation vector, identified 
as the new pion field. This construction also serves as a basis for the formulation of chiral 
perturbation theory (ChPT) [8], which is now considered as the realistic effective theory of 
low-energy hadron dynamics. In the context of the nonlinear realization of chiral symmetry 
a light cr-meson might be treated as a resonance in the TTTT-system [9,10]. Alternatively, for 
the linear case the a can either be described as a member of a 4-quark multiplet [2] or as 
quark- ant iquark resonance [11]. The different mechanisms for generating a light a do not 
necessarily exclude each other, but could in turn be related in a way which is not completely 
understood yet. 

Recently, the E791 Collaboration at Fermilab [12] and the BES Collaboration at 
BEPC [13] reported on evidence for a light and broad scalar resonance in nonleptonic cas- 
cade decays of heavy mesons. In the Fermilab experiment it was found that the a meson 
is rather important in the D meson decay D — Svr [12]. In a coherent amplitude analysis 
of the Stt Dalitz plot the scalar resonance is determined with 478^11 ^ MeV and total 
width 324:til ± 21 MeV. A fraction / = (46 ± 11)% of the decay mode D+ 7r+7r-7r+ is 
generated by the intermediate a-resonance channel. The measured branching ratio of the 
two-body decay air^ relative to the uncorrelated Stt decay of the D meson is then 

deduced as r(a7r+)/r(7r+7r+7r-) = 0.695 ± 0.135 ± 0.032 [1,12]. The BES experiment [13] 
concentrated on the nonleptonic decay J/^' — >• era; — > tttto;. The extracted values of the a 
mass and width are: = 390l|^ MeV and = 2821^^ MeV. 

Preliminary analyses of these two experiments were performed in Refs. [14-16]. In 
Ref. [14] the relevant coupling constants of the respective two-body decays D cttt and 
a 7171 were extracted from the data of the E791 experiment [12]. A direct calculation of 
the D — > (TTT amplitude was done in Ref. [15] in a constituent quark-meson model. Both 
analyses neglect the intrinsic momentum dependence of the D ^ aTT transition form factor 
and, in the case of Ref. [14], the final state interaction in the three-body decay D Stt. 
The two approaches [14,15] arrive at a disagreement between the analysis of the nonleptonic 
two- and three-body decays of the D meson. The extracted [14] or calculated [15] coupling 
constant goan is approximately twice as large as the one deduced from experimental data on 
the two-body decay D — > cttt [1]. In Ref. [16] the effective coupling constant gj/^auj was esti- 
mated using the perturbative QCD technique. Also, the role of the light cr as a elementary 
particle [17] and as a correlated two-pion state [18] was examined in B — > p7r decay. 

In the present paper we consider the two-body nonleptonic decays involving the light 
cr-meson with D — > crvr and B — > crvr. We work in the framework of a QCD motivated, 
relativistic quark model which implements a linear realization of chiral symmetry [19]- [21]. 
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In this context the formahsm also aUows to describe the pion as a composite particle. To 
solidify and extent our considerations, we therefore also present a comprehensive analysis of 
the electromagnetic form factors of tt and a associated with the transitions tt — > 7r7, tt — > 77 
and a — > 77. 

The specific scheme we work in can be viewed as an effective quantum field theory 
approach based on a Lagrangian of light and heavy hadrons (both mesons and baryons) 
interacting with their constituent quarks [19]- [21]. The coupling strength of a specific hadron 
to its constituent quarks is determined by the compositeness condition = [22,23], 
where Zh is the wave function renormalization constant of the hadron. The compositeness 
condition enables to relate theories with quark and hadron degrees of freedom to effective 
Lagrangian approaches formulated in terms of hadron variables only (as, for example, Chiral 
Perturbation Theory [8] and its covariant extension to the baryon sector [24]). 

Our strategy is as follows. We start with an effective interaction Lagrangian written 
down in terms of quark and hadron variables. Then, by using Feynman rules, the (S-matrix 
elements describing hadron-hadron interactions are given by a set of quark diagrams. The 
compositeness condition is sufficient to avoid double counting of quark and hadron degrees 
of freedom. The Lagrangian contains only a few model parameters: the masses of light 
and heavy quarks, and scale parameters which define the size of the distribution of the 
constituent quarks inside the hadron. This approach has been previously used to compute 
exclusive semileptonic, nonleptonic, strong and electromagnetic decays of light and heavy 
hadrons [19]- [21] employing the same set of model parameters. 

In such a way we consider the cr-meson as a quark-antiquark bound state of the light u 
and d flavors. We assume that the coupling strengths of both the pion and the sigma meson 
to the constituent quarks are identical in accordance with the linear realization of chiral 
symmetry. Based on this scheme the a meson mass is determined as — 385.4 MeV, 
which is in good agreement with the BES- result. Next we discuss in detail the gauging 
of the nonlocal meson-quark interaction and, for consistency, calculate the electromagnetic 
form factor of the pion. We then proceed with the 7r° — > 77 and a 77 processes. Here 
we obtain explicit expressions for the relevant form factors and evaluate the leading and 
next-to-leading order for large Euchdean photon virtualities. As a result we also obtain the 
two-photon decay width of the cr-meson. We compute the strong decay width for the process 
a Tin and flnally turn to the nonleptonic decays D (Ttt and B (Ttt. All our results 
are compared to available experimental data and to other theoretical studies. 

The layout of the paper is as follows. In Sec. II we begin by introducing the relativistic 
quark model which implements a linear realization of chiral symmetry. Sec. Ill is devoted to 
the derivation of the electromagnetic properties of the tt and a mesons. In Sec. IV we examine 
the strong and nonleptonic decay characteristics of the transitions a — > nn and D{B) an. 
Then in Sec. V we turn to a numerical analysis of the processes considered and discuss the 
quality of results in comparison with experiment and previous model calculations. Finally, 
we summarize our results in Sec. VI. 

II. THE MODEL 

We will consistently employ the relativistic constituent quark model [19]- [21] to compute 
a variety of observables related to the tt and a mesons. In the following we will present details 
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of the model which is essentially based on an effective interaction Lagrangian describing the 
coupling between hadrons and their constituent quarlcs. 

The coupling of a meson H{qiq2) to its constituent quarks qi and q2 is set up by the 
Lagrangian 

^nti^) ^ 9hH{x) JdxiJdx2FH{x,xi,X2)q2{x2)rHXHqi{xi) +h.c. (1) 

Here, Xh and Th arc Gcll-Mann and Dirac matrices which inclose the ffavor and spin quan- 
tum numbers of the meson field H{x). The function Fh is related to the scalar part of the 
Bethe-Salpeter amphtude and characterizes the finite size of the meson. To satisfy trans- 
lational invariance the function Fh has to fulfil the identity Fh{x + a,Xi + a,X2 + a) — 
Fh{x, xi,X2) for any 4-vector a. In the following we use a particular form for the vertex 
function 

Fh{x, Xi,X2) = S{x - W21X1 - Wi2X2)^Hiixi - ^2)^) (2) 

where is the correlation function of two constituent quarks with masses mi, 1712 and 
Wij = ruj/ {rrii + mj). 

The coupling constant qh in Eq. (1) is determined by the so-called compositeness condi- 
tion originally proposed in [22], and extensively used in [19]- [21], [23]. The compositeness 
condition requires that the renormalization constant of the elementary meson field H{x) is 
set to zero 

ZH^l-^fl'H{Mf,)^0 (3) 

where 11^ is the derivative of the meson mass operator. To clarify the physical meaning 

1/2 

of this condition, wc first want to remind that the renormalization constant is also 

interpreted as the matrix element between the physical and the corresponding bare state. 
For Zh = it then follows that the physical state docs not contain the bare one and is 
described as a bound state. The interaction Lagrangian of Eq. (1) and the corresponding free 
parts describe on the same level both the constituents (quarks) and the physical particles 
(hadrons) which are supposed to be the bound states of the constituents. As a result of 
the interaction, the physical particle is dressed, i.e. its mass and wave function arc to 
be renormalized. The condition Zh = also effectively excludes the constituent degrees of 
freedom from the physical space and thereby guarantees that double counting for the physical 
observable under consideration is avoided. Now the constituents exist in virtual states only. 
One of the corollaries of the compositeness condition is the absence of a direct interaction 
of the dressed charged particle with the electromagnetic field. Taking into account both the 
tree-level diagram and the diagrams with the self-energy insertions into the external legs 
(that is the tree-level diagram times Zh — 1) yields a common factor Zh which is equal to 
zero. We refer the interested reader to our previous papers [19]- [21], [23] where these points 
are discussed in great details. 

The meson mass operator appearing in Eq. (3) is described by the Feynman diagram in 
Fig. 1. In the case of pseudoscalar (Th = ij^) and scalar (Th = I) mesons, relevant for the 
present paper, we obtain the expression 
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2p2 o?p" J 47r2i 



rH'S'l(/(f + W21 ^)THS2{fi - Wi2 f) 



(4) 



where /c^) is the Fourier-transform of the correlation function ^ Hijy^x — ^-if) and 
is the quark propagator. We use free fermion propagators for the valence quarks 

m) = (5) 

with an effective constituent quark mass mj. As discussed in [19]- [21] we assume for the 
meson mass Mh that 

Mh <mi+ 1712 (6) 

in order to avoid the appearance of imaginary parts in the physical amplitudes. The cal- 
culational technique for determining the explicit expression of 11^ (p^) (4) is outlined in 
Appendix A. 

The interaction with the electromagnetic field is introduced by two ways. The free 
Lagrangians of quarks and hadrons are gauging in a standard manner by using minimal 
substitution: 

d^H^ {d^ ^ ieA^')H^, d^^q {d^ - ieA^')q, d^q {d" + ieA'')^, (7) 

where e is a proton charge. It gives us the first piece of the electromagnetic interaction 
Lagrangian written as 

CZ^'\x) = eqix) AQqi^) (8) 
+ ieA^{x) {H-{x)d^'H+{x) - H+{x)d''H-{x)) + e^Al{x)H-{x)H+{x). 

The gauging nonlocal Lagrangian in Eq. (1) proceeds in a way suggested in [26]. To 
guarantee local invariance of the strong interaction Lagrangian, in one multiplies each 
quark field q{xi) with the gauge field exponentional that gives 

^str+em(2)^^^ = ldx,JdX2FH (x, X^, X2)q2{x2) c''"^'^^''^^''^ (9) 



where 



I{xi,x,P)^ J dz^Ai^iz). (10) 



It is readily seen that the full Lagrangian is invariant under transformations 

qi{x) ^ e'^^i^^'^kiix), qi{x) qi{x)e-'''>if^''\ H{x) e'^"^^''^ H{x), 
A'^ix) A^'{x)+d>'f{x), 

where en — — e^^ . 
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Then the second term of the electromagnetic interaction Lagrangian arises, when 
expanding the gauging exponential up to a certain power of A^, relevant for the order 
of perturbation theory and for the process we consider. It seems that the results will be 
dependent on the path P which connects the end-points in the path integral in Eq (10). 
However, we need to know only derivatives of such integrals under calculations within the 
perturbativc series. Therefore, we use a formalism [19,25,26] which is based on the path- 
independent definition of the derivative of I{x, y, P): 

/j^n^^^A^^^'^'^) ^ _^^ira [I{x + dx,y,P')-I{x,y,P)] (11) 

where path P' is obtained from P when shifting the end-point x by dx. Use of the definition 
(11) leads to the key rule 

^ I{x,y,P)=A,{x) (12) 



which in turn states that the derivative of the path integral I{x,y,P) does not depend on 
the path P originally used in the definition. The non-minimal substitution (9) is therefore 
completely equivalent to the minimal prescription as evident from the identities (11) or (12). 
The method of deriving Feynman rules for a non-local coupling of hadrons to photons and 
quarks was already developed in Refs. [19,26] and will be discussed in the next section, where 
we apply the formalism to the physical processes considered here. 

For example, the piece of the Lagrangian in Eq. (9) in the first order over charge reads 

as 

^ 9hH{x) jdxijdx2 JdyE'^{x,xi,X2,y)Ai^{y)q2{x2)THXHqi{xi), (13) 

E^„{x,x^,X2,y) -J-^^J^J ^e^^^(--^)-^^(---)+^''(^-^)£:r(pi,P2,g), 

1 

EUpi,P2, q) = -eg,wniwi2q^ + 2p^) J dW^ {-tiwi2q + Pof - (1 - t)pl) 



1 

+ eq^W2iiw2iq'' - 2p^) j dWa [-t{w2iq - Pof - (1 - t)p\ 



PO = W12P1 + W21P2- (14) 

Transition matrix elements involving composite hadrons are specified in the model by the 
appropriate quark diagram. For example, the transition form factor which characterizes the 
hadronic transition Hi^^pi) — > if23(p2) + -^^12(^3) is determined from the Feynman integral 
corresponding to the diagram of Fig. 2: 

A^^^^^^^^(Pl,P2) = ^^H.3^...35^./^'^^^^^^(Pl,P2) 

l''-'''-'''{Pl,P2) = - l^^^„,,{-{k + W,,p,f)^H,,{-{k + W2sP2f) (15) 
X ^Hr,{-{k + Wi2Pl + W2lP2TMS2{l^+ 1^2)^' m,Si{l^+ ^i)r H.sSsil^)^ h,s] 
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where technical details concerning the derivation of integral (15) are indicated in Ap- 
pendix A. 

Finally we have to specify the vertex function (Eq- (4)), which characterizes the 
finite size of the hadrons. Any choice for is appropriate as long as it falls off sufficiently 
fast in the ultraviolet region of Euclidean space to render the Feynman diagrams ultraviolet 
finite. We employ a Gaussian for the vertex function ^^/(A;!;) = exp(— A;|;/A|^), where ks is 
an Euclidean momentum. The size parameters A|^ arc determined by fitting to experimental 
data, when available, or to lattice results for the leptonic decay constants fp where P — 
TT, D, B. The leptonic decay constant fp = Tp{Mp) is determined from [21] 



Mp')p' = ^ J^.^p{-k')tT\^0''Si{^ + W2i ^h'S^i^-w^^ f\ . (16) 

To reduce the set of values for A^/, we use a unified size parameter for hadrons with the 
same fiavor content. The best fit to the decay constants fp is obtained (see Table I) when 
the values of the constituent quark masses and the parameters A.h are choosen as follows 
[21] 



m,, 



u{d) 



0.235 0.333 



rrir 



rrih 



1.67 5.06 



GeV 



(17) 



and 



A 



7r(o-) 



Ak a 



1.6 



D 



A, 



2.25 GeV 



(18) 



In this paper we are aiming to explore the properties of pions and light cr-meson. We 
write down the relevant interaction Lagrangian in accordance with the two-flavor linear 
(7-model [3,4] 

^c^r^'^l^) -:^ldy Hy") + v/^) Wi^) + ^(^) r]q{x- y/2) . (19) 

Here the constraints on the bare couplings g = Q-k = Qa and the vertex functions $ = $^ = 
$0^ are imposed by chiral symmetry. We also require that the dressed couplings determined 
from the compositeness condition of Eq.(3) should be equal each other. This requirement 
allows us to determine the a meson mass by the equality 

%{Ml) = nUM^) , (20) 

where the physical pion mass M^^ is used as an input. 

With the hadron parameters fixed, from Eq. (20), we deduce a a mass of 

= 385.4 MeV. (21) 

The predicted value we obtain is close to the BES result [13] of — 390^3g MeV, which 
sets the lower scale for the range of mass values compiled in [1]. 
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III. ELECTROMAGNETIC PROPERTIES OF tt AND a 



In this section we proceed with the formahsm on the electromagnetic properties of tt and 
a mesons. For consistency we also include the tt meson in the discussion, where, in addition, 
the model predictions will serve to solidify the validity of the previously outlined approach. 
In the analysis we consider the following related processes: the electromagnetic transition 
of charged pions, that is vr"^ 7r'*'7, and the form factors characterizing the transitions 
7r° 77 and cr — > 77 for different kinematical regimes of the photons (real and virtual) All 
amplitudes considered are obtained in a manifestly gauge-invariant form. 

When gauging (9) the nonlocal strong interaction Lagrangian additional "contact" 
vertices are generated, which couple hadrons, quarks and photons, as already discussed in 
Sec. II. In particular, for the process tt"^ 7r^7 we need a vertex describing the coupling 
of a charged pion, a single photon and two quarks (Fig. 3a). The coupling of cr to a photon 
and two quarks (Fig. 3a) and to two photons and two quarks (Fig.3b) will contribute to the 
transition a — > 77. 

The full set of Feynman diagrams for the electromagnetic transition amplitudes consid- 
ered are summarized in the following: The transition vr^ 7r^7 is described by a triangle 
diagram (Fig. 4a) and the two additional contact diagrams of Figs. 4b and 4c. As mentioned 
in Sec.II, the diagram describing the direct coupling of pions to a photon is compensated 
by the counterterm (see Eq. (8). The process 7r° — > 77 obtains a contribution by the single 
diagram of Fig. 5a. The transition amplitude cr — >■ 77 is generated by the four diagrams of 
Fig. 5: a triangle diagram (Fig. 5a) and three contact diagrams (Figs.5b-d). In the follow- 
ing we denote diagrams containing a contact vertex with a single photon line (Figs. 4b, 4c, 
5b and 5c) as "bubble" -diagrams and the amphtude of Fig.5d as "tadpole" -diagram. The 
dominant contribution to the electromagnetic form factors arises from the leading triangle 
diagrams (Fig. 4a and 5a). We will demonstrate that diagrams containing contact vertices 
give negligible contribution, but in general they should be kept to guarantee electromagnetic 
gauge invariancc. Feynman rules for the evaluation of the nonlocal vertices of Figs. 3a and 
3b are derived in Appendix B. Next we will give the definition of the amplitudes for the 
processes tt^ — > 7r^7, 7r° — > 77 and cr —> 77 and discuss the properties of the relevant form 
factors. 



A. The electromagnetic form factor Ft^{Q'^) of the pion 

It is convenient to write down the vertex function for the transition 7r^(p) T^^ip') + 
7(g) in the form 

^'{P.P') = ^ P.(P^) - t,{p")] + A![(P,P') (22) 

from which the Ward-Takahashi identity follows immediately 

g^A'^(p,p')=S.(/)-S.(p'^). (23) 

For setting up the electromagnetic vertex function of Eq. (22) we use the pion mass operator 
E^(p2) with 
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47r2 



(24) 



^-k^)D{p) , = tr 



and A^(p, p') is the part which is orthogonal to the photon momentum: q^-^-j_{p,p') = 0. 
The exphcit expression for A^(p,p') results from the sum of the gauge- invariant parts of the 
triangle (A) (Fig.4a) and of the bubble (bub) diagrams (Figs. 4b and 4c): 

The separate contributions I'^^ (p, p') and /bub± (P^ P') given by 



i'LMp') 



$ - 



k + 



P 



$ - 



k + 



p 



tr[^'S{k + p')^^S{k+p)^'S{k)], 



IL^Ap^p')-K J ^*(-^') Jdt^'(^-k'-t[kq + ^]^ kri [d{p) - D{p')Y (26) 



77^ y AnH 



where ^'{z) = d^{z)/dz. Also, = 7/^ — gf/g^ and 77'* = P^ — {Pq/q^)q^ are orthogonal 
to the momentum transfer q with P = p 
In the limit g = we obtain 



^p^' 



dp'^ 



while, by definition, we also have 



A'^(p,p) = 2p'^F,(0) 



(27) 



(28) 



where F7r(0) is the pion charge form factor at the origin. Prom the comparison of Eqs. (27) 
and (28) it follows that the compositeness condition (3) is equivalent to the normalization 
of the pion form factor at the origin with Ft^IO) = 1. 

The full electromagnetic form factor F^^^Q'^) of the pion is defined by 



A^(p,p'; 



^pf^F^Q^) 



(29) 



where Q"^ — —q^ — —{p—p'y is the Euchdean momentum transfer squared. The elec- 
tromagnetic radius of the pion, related to the slope of Ft^{Q^) at the origin, is then given 

by 



< r >= 



-6 



dF^Q^) 



(30) 
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B. Transitions tt 77 and cr 77 

The matrix elements for the transitions 7r° — > 77 and cr — > 77 can be written in general 
in the manifestly gauge-invariant form 

52) = e' a^'' F^,,{p', qf, qj) , (31) 

Mr-.,,iqi, ?2) = I 6'''^ F,^^(p^ ql qj) + c^^ G,,,{p\ ql qj) } . (32) 

The tensors a'^", h^^ and c^" refer to the Lorentz structures: 

= e^-/^ (g,), {q^)p , 

ft"'^ = ^^'^ (?i?2) - grg2^ , (33) 
c^'^g^Ulql + ?^2'(gig2) - grgr?2 - , 

while gi and 52 are the photon four-momenta and p — qi + q2is the meson momentum. Here, 
Fir-y-yiP^ 1 Ql, Q2) is the tt" — > 77 form factor. The transition cr — > 77 is characterized by two 
form factors: Fo-^^(p^, qf, g|) and Ga-yy{p'^, qf, q^)- Usually only one of the form factors, that 
discussed in the literature. However, when both photons are off-shell, constraints 
set by gauge invariance result in the two terms of the matrix element M^^^^, proportional 
to the Lorentz structures b^^ and c'*''. If at least one of the photons is real, then only the 
first form factor F^-y-y gives a non- vanishing contribution to the invariant matrix element. 
The decay width of the transition if — > 77 with if = tt or o" is given by 

Th^^ = ^a'Mjjgl^^ (34) 

where gn-yy = Fn-yyiM'^, 0, 0) is the 7777 coupling constant and a = e^/(47r). 

For the evaluation of the ii — 77 form factors we pursue the following strategy. 
First we consider the calculation of the form factors Ft^^-^{p^ , q^, q^) , F^^^ij? ,qf,q2) and 
Ga'yyip^, ifiOi) in so-called local limit. In this case the dimensional parameter A^, appear- 
ing in the 77(0") meson vertex function $(— A;^), is taken to infinity or, which is equivalent, 
$(— /c^) — > 1. Wc obtain analytical expressions for the if — 77 form factors and show that 
for the case of 7r° 77* the form factor has the incorrect behavior for large photon virtuali- 
ties in the Euclidean region. Then we repeat the derivation using the nonlocal approach and 
show that the correct asymptotics for the 7r77* form factor can be reproduced when using 
the mesonic quark- antiquark wave function. The final numerical analysis for the complete 
form factors will be done in a separate section. 

1. Local 7r°77 and (777 diagrams 
The Feynman integral corresponding to the local (L) triangle diagram is written as 

IZ:^{<11,<12) = /^tr{7^5(^+ ^i)rM(^- ^2)7^5(^)} (35) 

where F^r = i7^ for vr and F^- = i for a. For simplicity we drop a common factor of — 3/47r^. 
The complete integral (35) is free of the logarithmic ultraviolet (UV) divergence. Separate 
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contributions to Eq. (35) with the UV divergences are treated using dimensional regular- 
ization [29]. In this case we use the definition for 7^ in D dimensions (see details in [30]) to 
guarantee the fulfilment of two important properties: (7^)^ = / and tr(7^7'*7^) = 0. For the 
pseudoscalar triangle the UV divergency is completely gone due to the contraction of the 
divergent integral with the trace tr(7^7'*7^), which by definition is zero in any dimension. 
In the case of the 0-77 diagram the use of the gauge- invariant regularization (e.g., dimen- 
sional regularization) guarantees the cancellation of the UV parts and generates the correct 
finite part. The finite part of the 1777 diagram is the same for both types of gauge-invariant 
regularizations. The master formulas used in dimensional regularization are 

Ub'Yl') I d-'k , = (for meson) 



/ 



[w? — k'^ — Dq 

4A;^r + g^'-(m^-fc^-Do) 

a A; 1-^ — 713 — = (for a meson) 

[m^ — — DqY 



where Dq — aia2P^ + aia^qf + Q!2Q!3?2 with p = qi + q2 arises when using the a- 
parametrization. 

The integral Eq.(35) is gauge invariant 

^Z'Li(li^ q2) ■ {ql)^l = Ia,l{(1u 0.2) ■ (^2)1. = 
as can be easily deduced from the Ward identity 

and simple algebra. 

The integral of Eq.(35) is particularly simple for the case of the pion with 

I'Z'^iQu ?2) = m a^^ Flip', ll qI) (36) 
where the form factor Fl is given by the two-dimensional integral /l: 

Fl{p\qlql)^I^{p\qlql)^ [ d'aS (l -^.a] ■ (37) 



=1 



The double integral Jl can be further reduced to a single one using the 't Hooft-Veltman 
technique [31] and can be finally expressed by a combination of Spense functions. Here, 
however, we will stay with the integral representation. We have 

Hp ,q,,q,)=Jdxl ^y,.^^ 2 2 (38) 



qi - oiQp 

u 

h\[l — x{l — x) ql/m^] ^ h\[l — x{l — x) ql/rTn?] 



where 



AV2 .x-ql-ql q;o/(1 - ckq) AV2 .x + ql{l- q;o)/q;o + ql 



A = X{p\ ql ql) ^p' + qt + qt - ^pVi - 2pV, - 2qlql (39) 
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is a kinematical triangle function and = {jp + qf — ^2 + A^/^)/2p^. All three thresholds 
can be readily seen from this representation. In particular, when both photons are on their 
mass-shell (g^ = g| = 0) we get 



^ , o , 1 /■ , ln[l — a;(l — a;)p^/m^ 



(40) 



dv In 



J — — i 

Am? 



1 



1 + v^r^ \ 

1 — \/l — V J m? — V ■ / A — ie 



— le p^ 



Aw? 



where p{k'^) = In { (l + ^l-Am?/K^) / (l - ^l-Am?/^'^) } is the spectral function. 

Another interesting limiting case is when one of the photons has large Euclidean mo- 
mentum, for instance, q\ — —Q^, ?2 — p^ — 0. In the limit — > oo the integral 
Eq.(38) reduces to 



/L(o,-g^o) 



dx 



ln[l + x{l - x)Q'^/w?] \n^{Q'^/m 



x{l — x) 



2g2 



Note that such an asymptotics is in contradiction to the QCD-prediction for the 7r°7*7 form 
factor of 1/Q2 [32]. 

In the following we turn to the evaluation of the CT77 diagram in the local limit using 
dimensional regularization. The defining integral of Eq. (35) can be written as 



(41) 



with 



FZip\qlql) 



Gl{p\qlql) 



{I-A,L ■ b) - {ILl ■ c) jb-c) 

(/i,L-c)&^-(/X,L-&) {h-c) 
62c2 - (6 • C)2 



(42) 



and where the dot-product refers to the contraction of the Lorentz indices. 
Using the a-parametrization we have 



jd^a5{l-j2 



OLi 



m? — Dq ' 



(43) 



m 



Glip^qlqi)^ 

HI (I2 



^3 ./^i ^^ \ ai(l-2ai)g2 + a2(i_2a2)gi 
d ad \1 — ai I 



— Dq 
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Again, one integration in Eq. (43) can be performed analytically (the resulting expres- 
sions (106) and (107) are indicated in Appendix C). For the case qf — q2 — we get 



1-V 



w? — V ■ ' 

(44) 



'^^(^■''•°)-4j^/*{i('^''""(l 



— ^/l — V J j — V ■ p'^/4: 



2. The nonlocal 7r°77 and C777 diagrams 

For a nontrivial, that is nonlocal, meson-quark vertex the evaluation of the relevant 
Feynman diagrams is based on a method outlined in [19]. This technique was originally 
developed to derive a representation for the 7r'^77 diagram (Fig. 5a) with a dressed pion- 
quark vertex, which is described by an arbitrary function $(— /c^) decreasing rapidly in the 
Euclidean region. It also preserves translational and gauge invariance in a manifest way. 

In our model the 71^77 form factor is given by 

F^^^(p^ ql, ql) = ^J^^ ^ -^NL(p^ ql, ql) (45) 

where Jnl is the nonlocal (NL) quark-loop integral: 

d'^k l>(-A;2) 



-^nl(/, ql,qP 



At:H [m? -{k+ p/2)2] [m^ - {k - [m^ - {k + q/2f] 

(46) 



00 2 / 3 \ 



with q = q2 — qi- 

The argument z appearing in the vertex function is written in the form 



f „ t 



D + A, 

1+t 1+t ' 

D — m — aia2P — aia^ q^ — a2a3 q2 = m — Dq , 
2 

A = m - J + Y (P - ?i - ?2) ■ 

We choose to define the functional dependence of the vertex function as $(— /c^), hence after 
transition to the Euchdean region (A;° — > iki or k'^ — > —k^) the argument changes sign. 
This convention allows a consistent choice of the functional form for $(fc|;) in the Euclidean 
region, where we finally perform the integrations. For example, particular forms of the 
vertex function are ^{k^) = exp(— /clj/A^) (Gaussian vertex function) or $(A;|;) = 1/(1 + 
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k'^/h?)^ (pole vertex function), where 1/A characterizes the size of the meson. Obviously, 
the limiting local case /nl II can be recovered for A — oo. This limit also serves as a 
check for the numerical calculations used in the final step of the evaluation. The following 
analytical results are obtained for an arbitrary form of l>(— fc^). Only when turning to the 
final numerical calculation, a specific form of ^{k%) (Gaussian vertex function) will be used. 

Again, one integration in Eq. (46) is performed by linearization of the argument z with 
respect to one of the a parameters as it was done in the local case. The resulting expressions 
are indicated in Appendix D, as given by Eqs. (109) and (110). For the study of the 
dependence of the 7777 form factor on the photon virtualities in the Euclidean region it is 
useful to introduce the variables: qf — —(1 +0;) and g| = —(1 — a;) where Q"^ and 
uj are the total virtuality and the asymmetry, respectively^ . With this convention we define 
the pion form factor as 

F^,.,.{Q',u;) ^ F.^^(m^,-(1 + (^)^,-(1-(^)^). (47) 

The particular choice cu — ±1 corresponds to the physically interesting case where one of 
the photons is virtual and the other is on its mass-shell. The corresponding transition form 
factor F^^^*(Q^) = F^^^(M^, —Q^, 0) has recently been measured by the CLEO Collabora- 
tion [33] for momentum transfers in the range from 1.5 to 9 GcV^. A discussion of previous 
experiments can also be found in Ref. [33]. The CLEO data confirmed the predictions of 
perturbative QCD (pQCD) [32] for the asymptotic behavior of F^-y-y*((5^) at large Q'^: 

K^AQ^) = ^ + 0(i,) (48) 

where = f^/\^. 

A detailed analysis of the different QCD approaches to the F^^^. (Q^) form factor has 
recently been done in [34]. Unfortunately, a straightforward calculation of F^^^*((5^) in the 
context of QCD is either not fully possible, because the operator product expansion (OPE) 
fails at small Q^, or contains unknown parameters like the twist- four scale parameter 5^ re- 
lated to the gluon condensate [35,36]. Therefore, an analysis of -Ftt^^* in our QCD-motivated 
approach seems to be quite reasonable. Information about the i^f — 77 form factors for 
non-trivial photon virtualities is relevant to forthcoming experiments on production of pseu- 
doscalar mesons in 7*7* collisions. Also, there is the project to search for a hght scalar 
meson in very peripheral heavy ion collisions (for a recent review see Ref. [37]). 

For the analysis of F^^Y'r* perform a systematic expansion in powers oil/Q^ including 
leading order (LO) and next-to-leading order (NLO) terms with 

....(«^.)^...{^.^.o(±)}. 

The expansion coefficients J^^{u) and J^^'~'{u) are given in our approach by 



^Another possibility is to work in terms of the average virtuality. In this case the variable 
should be rescaled as — > 2Q^ in all further formulas. 
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3 i?LO(0) 



J. 



NLO, 



NLO/ 



2 R 



3 i?LO(o) 



(50) 



The structure integrals R}f{uj) and i?jj^'"°(a;) are derived in Appendix D applying the chiral 
limit with Ml — and are written as 



R. 



LO 



dt 



j^^{mH) Ln_(i,u;), R. 



,NLO 



(o;) = - 4 y^dt l>(m^i) D(i, u) , (51) 





where 



Ln±{t, Lu) = In 



1 + t(l +a;) 
1 + t 



± In 



l + t{l-u) 

T+t 



[l + t{l + u)] [l + t{l-u)] 



In the derivation of Eqs. (49) and (50) we use the identity 



87^2^/2 " 



LO 



(0), 



(52) 



relating and the meson-quark coupling constant g. From the expressions of Eq. (51) it 
can be readily seen that the expansion coefficients contain only even powers of uj, which is 
a consequence of Bose-Einstein symmetry and charge conjugation invariance. 

At this stage of the development we indicate a first comparison of our results for the ex- 
pansion coefficients calculated at a; = 1 and a; = to the ones predicted by pQCD [32,35,36]. 
A summary of previous model results for J^*^(l) and J^*^(0) can be found in Ref. [38]. First, 
we consider the experimentally accessible case w = ±1, where one of the photons is real and 
and the other one is virtual. Our numerical result for J^^(l) = 0.998 ~ 1 is in very good 
agreement with the pQCD prediction of J^°(l) = 1 [32]. The deviation of our value for 
J^°(l) from 1 is explained by the use of simplest form of quark propagator and meson-quark 
vertex form factor. 

In the relativistic quark model the NLO coefficient J^'"*^(l) is proportional to the con- 
stituent quark mass squared and results in J^^*-'(l) = —0.37 GeV^. The pQCD approach 
predicts [34] 



80 2 
-27^ ' 



(53) 



where 5^ is the twist-four scale parameter related to the gluon condensate [35,36]. Originally, 
this quantity was estimated to have a value of 5^ = (0.2 ± 0.02) GeV^ using the QCD sum 
rule approach [36]. A recent evaluation [34] results in: 5^ = (0.19 ± 0.02) GeV^. By taking 
into account the value of Ref. [36] the pQCD prediction is 



jNLO(l) = -(0.59 ±0.06) GeV" 



(54) 



Now we turn to the case uj = 0. We exactly reproduce the LO coefficient predicted by 
pQCD [36]: J^°(0) = 2/3. For the NLO coefficient our result is 
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JNLO(O) = -Im^ = -0.074 at m = 235 MeV . (55) 
The result of pQCD is again proportional to 5^ [36] : 

JNLO(O) ^ - 2±5^ ^ -(0.24 ±0.02) GeV^ (56) 

Prom this comparison we conclude that our prediction for the LO expansion coefficients 
J^°(l) and J^°(0) are in perfect agreement with the expectation of pQCD. The NLO 

coefficients are rather sensitive to the choice of the constituent quark mass. In particular, 
when using m = 395 MeV in the calculation of J^^*-'(l) and m = 420 MeV in the case of 
jNLO^g-j ^Yie central values of the NLO coefficients predicted by pQCD can be fitted. 

For the process cr ^ 77 we have to evaluate, in addition to the triangle diagram (Fig. 5a), 
the bubble (Figs. 5b and 5c) and tadpole diagrams (Fig.5d), which arise from gauging the 
nonlocal aqq interaction Lagrangian. Each particular diagram is not gauge invariant by itself, 
but the total sum fulfills the gauge invariance requirement. To simplify the calculation we 
split the contribution of each diagram into a part which is gauge invariant and one which 
is not. This separation can be achieved in the following manner. For the 7-matrices and 
vectors with open Lorentz indices /i and u we use the representation: 

(57) 

7 - 7±;,2 + 92 -2 > ^ = ^±;«2 + ^2 -2- > 

92 92 

such that 7^.qj (gi)^ = ^.^^ (gi)^ = and 7^.^^ (^2)1. = w^jg^ (92)1/ = 0. Expressions for 
diagrams containing only ±-values are gauge invariant separately. It is easy to show that 
the remaining terms, which are not gauge invariant, cancel each other in total. 

The gauge-invariant parts of the nonlocal (777-triangle (A), bubble (bub) and tadpole 
(tad) diagrams are given by 



-'A, 



= Fl^ (/, 9?, 92) + (/, 9?, 92) c^"" , (58) 

Cbj9i,92) = - / ^Jdr^'{-x{0,q,)) {kL^^^tr [fl,,,S{l^+ i,/2)S{}t- i,/2)\ } 



+ (?i ^ g2, - ^) = Kubx ip\ qI <&) b"'' + C^ubx qI <&) c'"" , (59) 

ATdj9i,92)=/^tr[5(^)] 'jdr [- ^ ■ {4.'(-x(0,p)) + 4.'(-a;(0, g))} 
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1 

+ r Jdl ( k + 





92 



^ / J_;gi 



^1;.. \ ^"i-<qi, (12)) + 52)) 

+ (91 ^ 92, ^ ^) - 92) ^'^^ + GLi, (P', ql ql) c^" , (60) 



where 



2^(91,92) = k"^ + kT{lqi + 52) + -^{Iql + 2lqiq2 + ql) . 



Applying the a-parametrization for the gauge-invariant part of the triangle diagram gives 
Fl^--^Jdt-^^^—^Jd^as(^l-p^aij {-l>'(^)} [t (1 - 4 ai as) + 2 as] , (61) 



G 



A, 



X 



m 

qlql 



m 



2 / 3 \ 

^^(Y^/^'«^(l-E«^J {-^'(^)} (62) 
(1 + 2 ai t)(l + 2 ao ^) 9 t (ai — 0:2) x 9 / ^ v 9\ 1 



Again, using the t'Hooft-Veltman technique one integration in Eqs. (61) and (62) can be 
performed. The resulting expressions are indicated in Eqs. (117) and (118) of Appendix E. 
The analytical results for the form factors F{G)^^y^^ and -F(G)5^adx given in Appendix F 
(Eqs. (125)-(128)). 

In the local limit, that is A —> 00, we obtain 



Fa, 



L 5 



m 



^bubx 



2 9192 

We therefore recover the local (T77 form factors: 



m 

2 2 1 
9i 92 



0, 



m 



2 9? 92 



-f'NL - i'A, + -f'bubx + -^tadx " 
^NL — ^Ax + ^bubx + ^tadx 



Fl, 



which were already indicated in Eq. (43). 

The form factors F,^^^ and Go-y-y are expressed in terms of the functions introduced above 

as 



-^a77(P^ 91,92) = 
Gaj-yip'^, 9l, 92) = 



59 



^9 



^Xx (P'' 5?' ^2) + Kubx (P'' ?i ' 52) + ^tW (^^'' 5i ' 52) 
Gax(/,9?,92') + Gl^^Jp^qlql) + Gl^^{p\qUl) 



, (63) 
■ (64) 



In analogy to the case of the pion we define the (777 form factors using the and uj variables 
accordingly: 
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F<.,.,.(g^a;) ^F.,,(M,^-(l + a;)^,-(l-o;)^) , (65) 
G,^.^.(g^a;) = G.^^(M,^ -(1 -(1 - o;)^) . (66) 

The final numerical analysis indicates that the additional nonlocal diagrams (bubble and 
tadpole) are significantly suppressed. For simplicity below we discuss the power expansion 
for Ffj'y*^*{Q'^ ,uj) and Gcr7*7*((3^,'-^) in the limit = 0, where only the gauge-invariant 
part of the triangle diagram is included. In the numerical analysis (Sec. V) we will take into 
account all diagrams and use the value of M^r = 385.4 MeV as predicted by our approach. 
For the power expansion of the (777 form factor we pull out, with the help of Eq. (52), a 
common scaling factor 2F^: 

F iO- u) - 2F 1'^'°^"^) + ^^'^'"(^) , o( ^ 
F^,.,.{Q ,cu) - 2F^ + + 0(^— 

(67) 



The expansion coefficients J^'~*(a;), J^^^(a;), K^^{u!) and K^^^{u!) are derived in the limit 
indicated above as 



9 i?LO(o) ' 



(68) 



where 



RTi^)--^ /rft$(m^t){^[l + i^(t,a;)] - ^ W(t,a;) 


^ ' 

Final numerical values for the expansion coefficients are: 

J^^il) = = 0.54, J,LO(0) = -^^^ = 0.23, 

J^^O(l) = 1.18 GeV' , ^^"^^(0) = 0.87 GeV' , (70) 

K^^O{l) ^ 4.61 GeV' , i^^'^°(0) = -0.09 GeV' . 
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IV. THE SIGMA MESON IN STRONG AND WEAK DECAYS 



In this section we derive the matrix elements describing the strong decay cr — > tttt and 
the nonleptonic decays D — > cttt and B — > air. 

A. Strong decay a tttt 

The strong form factor Gaivnip^, ^i) related to the transition cr ^ tttt is defined in our 
approach as 

>(-hf]l 

x.p.(.. 1).} 

The a momentum is given by p, pion momenta are labeled by qi and q2 with p = qi + q2 
and q — q2 — qi- The matrix element of Eq. (71) was already introduced in its generic 
form by Eq. (15), when we discussed the model features. The corresponding evaluation is 
summarized in Appendix A. 

Since the strong modes a — > 7t~^tt~ and a — > 7r°7r° dominate the a decays, the total width 
Fa is given by 

~ r((7 ^ TT+TT") + r((7 ^ TtV) . (72) 

With the coupling constant Qa-Kir = Ga-K-n-iM^, M^, M^) the total decay width can then 
be expressed as 

r. = ^ r(<7 ^ TT+TT-) = AV^(M,^ M^, M'^), (73) 

where A1/2(M2, M^, M^) was aheady defined in Eq. (39). 

In our estimate for the GainT form factor and the c-mcson width final state interaction 
is neglected. An accurate analysis of the strong decay properties of the a meson should also 
include these effects. We just refer to two Refs. [9,10] where the cr-meson was generated 
dynamically in the iso-scalar S-wave of tttt scattering in the context of chiral perturbation 
theory. The pole identified with the light a meson occurs ai E = Ma- — (i/2)rcr with 
Ma ^ 500 MeV and ^ 600 MeV [9] and 400 MeV [10]. Judging from this work, an 
inclusion of final state interaction will possibly lead to a further increase of as compared 
to the value derived in our model. 

B. Weak decays D ^ an and B ^ an 

Now we turn to the discussion of the nonleptonic two-body transitions H ^ an with 
H = D or B. The effective interaction Lagrangian relevant for the nonleptonic two-body 
decays D ^ an and B ^ an is given by [39-42]: 
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■•nl 
■'int 



^_Gf( 



a\0\ 



+ H.c. 



(74) 



where is the Fermi couphng constant and Vqqi are the matrix elements of the Cabibbo- 



Kabayashi-Maskawa quark-mixing matrix iy*^ Vud 
four- quark operators O*^ are defined as: 



0.217 and 



iCi)y-A{<ijd,^y^Ai 



01 = {diCi)v-A{ujdj)v-A, O2 = 

0\ = {Uihi)v-A{d3Uj)v-A, O2 = {dibi)v-A{UjUj)v-A, 



0.0036 ) [1]. The 

(75) 
(76) 



where i,j are the color indices and label V — A is a, short-hand notation for the 7'*(1 — 7^) 
Dirac structure. The couplings af are the combination of the Wilson coefficients including 
both the factorizable and the nonfactorizable effects. We use the following values for the 
effective couplings af: = 1.274 [41] and a\ = 1.038 [42]. 
The corresponding decay width r(i? — > air) is given by 



an] 



\'l\Ml,MlMll 



(77) 



where gHa-w is the effective weak coupling constant. The constant Qhu-k is equivalent to 
the expectation value of the effective Hamiltonian H-eff as derived from the nonleptonic 
Lagrangian (74) 



(78) 



where F^'^ is the weak form factor describing the H ^ a transition [14,15]. 

In Refs. [15,17] the authors included two contributions to the form factor F^'^ : the 
direct diagram of Fig.6a and the resonance or polar diagram of Fig. 6b with an intermediate 
axial meson H{1'^). The contribution of the resonance diagram is sizeable and close to the 
contribution of the direct diagram. As result they overestimated the experimental result for 
r{D"^ (T7r+). Therefore, we restrict ourselves to the consideration of the direct diagram 
(Fig. 6a) only. We neglect also the suppressed "annihilation" diagram, that is the D W 
a + TT transition, since it involves form factors of light mesons at high momentum transfer 
(g^ = M|)) [40]. Note, that the diagram generated by operator O2 is vanishing because the 
corresponding matrix element is proportional to < 0\{qq)Y-A\o' >= 0. 

The contribution of the direct diagram (Fig. 6a) to the form factor F^'^ is given by 



47r2V2(M^-M2) 



(79) 



where Ina is the structure integral: 



k + 



P 



X tr 



(80) 



with ujQq = tTiq/ {mQ + mq). Again, the analytical evaluation of the structure integral, before 
the final numerical calculation is applied, is indicated in Appendix A. 



20 



V. NUMERICAL ANALYSIS 



A. Electromagnetic form factors 

Preliminary model results for the pion charge Ft^{Q'^) and form factors were 

already presented in Ref. [19]. Here we extend our formalism to the case of the (777 form 
factors and we also perform a comprehensive analysis of the if — 77 form factors. 

Recently, new and more accurate experimental results for the charged pion electromag- 
netic form factor Ft^{Q'^) were obtained by the Jefferson Lab F^^ Collaboration [43]. These 
data for the momentum transfer region of = 0.6 — 1.6 GeV^ were extracted from an 
analysis of electro-production of pions on the nuclcon. The new results for the pion form 
factor lie somewhat higher than the older Cornell data points [44], but are consistent with 
a monopole parametrization fitted to elastic data at very low and with the scaling 
law [45]. In Fig. 7 we show our results for -^^(Q^) in the region up to = 4 GeV^. For 
comparison we indicate the experimental data, recent results (F^r Collaboration [43]) and 
previous ones (DESY [46] and CERN NA7 [47] Collaboration). We also indicate the predic- 
tions of other theoretical approaches: QCD sum rules [48] , light-cone quark model [49] , N JL 
model with a separable qq interaction [50] and QCD modeling approach based on solutions 
of the Dyson-Schwinger equations [51]. Our model predictions provide a rather good de- 
scription of the available data and are very close to the QCD sum rule result [48] including 
0{as) correction. They are also close to the results of the QCD motivated approach [51], 
which is based on a similar physical picture. 

In contrast to the work of Ref. [51] we use, roughly speaking, simple phenomenological 
prescriptions for the quark propagator (a free propagator with an effective quark mass 
instead of a confined one) and for the meson correlation function. We also get a reasonable 
description of the pion charge radius with = 0.65 fm, as obtained in our model. Our 
result should be compared to the present world average data of r^r = (0.672 ±0.008) fm from 
PDG2002 [1], to the recent experimental result = (0.65 ± 0.05 ± 0.06) fm of the SELEX 
Collaboration [52] and to the prediction of Ref. [51] with r^r = 0.67 fm. 

Next we discuss the numerical results for the 7r° — > 77 transition form factor. First 
we consider our results for the form factor F^^^. (Q^), which are given in Fig. 8. The data 
points are taken from [53] (CELLO) and [33] (CLEO). Other theoretical calculations include 
the hard scattering approach (HSA) [54], QCD sum rules [55] and pcrturbative light-cone 
QCD [56]. Our curve for Q^F^*^t^{Q^) is in good agreement with the data and approaches 
the Brodsky- Lepage limit for large Q"^. With the usual definition for the range of a form 
factor, that is 

with H^n,a, (81) 

g' = o 

we obtain for the radius of F^^^^* (Q^) : 

< rl^ >= 0.44 fm' . (82) 

Our result confirms the monopole- type approximation of the CLEO data [33] and is very close 
to the CELLO measurement [53] of < >= 0.42±0.04 fm'. Again, our model prediction is 
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close to the result of a similar theoretical approach with < r^^ >= 0.39 ±0.04 fm^ [38]. Fig. 9 
contains our results for the F^^.^. (Q^, cj) form factor for different values of the asymmetry 
parameter u = 1, 3/4, 1/2 and 0. An increase in uj leads to a rise of F^^.^. (Q^, ti;) for large 
Q^. Finally, for the coupling constant g-i^-yy and the decay width F^^^ our approach gives: 

(83) 



fl'7r77 



0.263 GeV" 



and 



F 



7r77 



7.15 eV 



which is close to the data [1] of g^^^^ = 0.273 GeV^"*^ and Tj^^^ = 7.7 ± 0.5 ± 0.5 eV. 

Now we turn to the discussion of the dependence of F^^*^*. First we consider the 
limiting case a; = 1, where one of the photons is virtual and the other one is real. The 
corresponding form factor F^^^*{Q'^) is plotted in Fig. 10. In the numerical calculation all 
three types of diagrams (A, bubble and tadpole) are included. As was already stated before, 
the A diagram gives the dominant contribution, whereas bubble and tadpole diagrams are 

form factor as 



significantly suppressed. Using Eq. (81) we determine the slope of the F^r-yy* 

< ri, >= 0.40 fm^ . 



IT7 



(84) 



For completeness we also present the results for Fcr^*j* for different values of the asymmetry 
parameter a; = 1, 3/4, 1/2 and in Fig. 11. 

The coupling constant g^-y-y and, according to Eq. (34), the decay width F^-^^ are given 
in our approach as: 



g^y^ = 0.330 GeV" 



and 



IT77 



0.26 KeV. 



(85) 



Again, the dominant contribution to g^^^^^ arises from the triangle diagram. The separate 
contributions of the bubble (bub) and tadpole (tad) diagrams to the coupling constant are 



g^^^ = - 0.4 X 10-^ GeV-^ and g^^^ = - 0.2 x 10"^ GeV"^ 

bub tad 



(86) 



A variation of the a meson mass in the region of < M^- < 2m = 0.470 MeV does not 
have much influence on the value for g^^^ with g„^^ ~ 0.31 ±0.02 GeV~^. In the soft meson 
mass limit — = we approximately reproduce the low-energy theorem: 

9_ 1 

To 



fl'7r77 — 



47r2F. 



(87) 



0.263 GeV-^ and g^^^ 
= 134.98 MeV and 



The constants gH-yy calculated for zero meson mass values, 5'7r77 ~ 
0.292 GeV~^, are rather close to the quantities predicted at M^- 
385.4 MeV. 

Taking the recent result for the two-photon decay width of a or /o(400 — 1200) of 
F(/o(400 - 1200) ^ 77) = 3.8 ± 1.5 KeV [1,57] at face value, our direct resuh of Eq. 
(85) is off by about an order of magnitude. However, predictions for the decay width de- 
pend rather sensitively, that is to the third power, on the value of the scalar meson mass 
as evident from Eq. (34). Using our canonical value of g^^^ — 0.330 GeV~^, which, as 
discussed, is fairly independent of the mass value, and varying the a mass we obtain: 

(88) 





GeV 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1 


1.1 


1.2 


F((7 ^ 77) , 


KeV 


0.3 


0.6 


1 


1.6 


2.3 


3.3 


4.6 


6.1 


7.9 



The range of prediction can be summarized as T[a — > 77) = 4.1 ± 3.8 KeV, which now is in 
qualitative agreement with the experimental result. 
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B. Strong decay a ^ tttt 



Based on the definition given in Sec. IV A, our numerical result for the coupling constant 
of the strong decay cr — > tttt is: 

^,,, = 1.8GeV. (89) 

The value we obtain is close to the one extracted from the BES experiment [16] with Qo-tt-k — 
2-0lo;i9 GeV and to the prediction of the hnear a model: 

W = ^ = l-6 GeV. (90) 

In Eq. (90) we used the values M^. = 385.4 MeV and = 92.7 MeV as determined in our 
approach. Using Eq. (73) for the decay width of the a meson we get r^. = 173 MeV. This 
value is smaller than the ones reported by the E791 (324t4o±21 MeV) and the BES {2S2tll 
McV) Collaborations. Again, a variation of the a meson mass in the region < M^r < 2m 
leads to a range of predictions for g^^^i^M"^) = 1.96 ± 0.73. The central value of gcnrniM^) 
corresponds to M„ ~ 350 MeV, the upper limit to = 0, whereas the minimal value 
is obtained for a mass value near M^. ~ 2m. Substituting the results for Qa-KniMl) as a 
function M„ into Eq.(73) we estimate the variation of the a meson width as Fg. = 206 
MeV. The lower limit corresponds obviously to the threshold — ^Mt^ and the largest 
value to Ma 330 MeV. 



C. cr meson in weak decays 



Finally we discuss the results for the nonleptonic decays D ^ an and B an. We first 
give our predictions for the weak form factors F^'^{M^) and Fq'^{M^) of Eq. (79) evaluated 
at the physical point — M^: 

F(f"(M^) = 0.298 and i^o^"(M^) = 0.141 . (91) 

With our results for F^"'{M^) and Fq'^{M^) we get for the effective weak couphng constants 
(Eq. (78)) 

QDan = 298 eV and gean = 15.8 eV (92) 

and finally for the decay widths, originally defined in Eq. (77), 

r(D+ (77r+) = 0.90 X 10-^^ MeV and F(S+ ^ (J7r+) = 0.94 x 10"^^ MeV . (93) 

Our prediction for the decay width F(D — > air) is in agreement with the lower value of the 
E791 result [1,12] of: 

r{D+ (J7r+) = (1.32 ± 0.31) x 10"^^ MeV . (94) 

We also compare our results to previous theoretical calculations done in Refs. [14,15,17]. A 
value for F^'^{M^) was estimated in Ref. [14] using the D ^ an ^ Sn data [12] without 
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properly taking into account the rescattering effects. The result is F^'^(M^) = 0.79 ±0.15, 
which is twice as large as our prediction and also results in a value for the width r{D^ — > 
cr7r+) larger than the measured one. Obviously, the value for F^'^{M^) should be directly 
extracted from the two-body transition D — > cttt. 

In Ref. [15] the form factor F^'^{M^) was calculated (including direct and polar contri- 
butions) in a quark model similar to ours. The contribution of the direct diagram 0.30 ±0.02 
is close to our result but the contribution of the polar diagram is very large 0.22to;oi. With 
the final result of r{D~^ — > cttt"'") = 2.3 x 10~^^ MeV this approach overestimates the exper- 
imental result. This work was extended in Ref. [17] to determine the form factor Fq^'^(M^). 
The result they obtain is F^''(M^) ~ 0.45 ± 0.15 which, again, is about three times larger 
than our prediction. 

VI. SUMMARY 

In summary, we have applied the rclativistic constituent quark model to investigate 
electromagnetic tt and a meson form factors and a variety of strong and weak decay char- 
acteristics involving the a. We start with an effective quark-meson interaction Lagrangian, 
which is based on a linear realization of chiral symmetry. Then we write down the matrix 
elements describing the meson interactions in terms of a set of quark diagrams. All model 
parameters were previously determined [19]- [21] in a fit to experimental observables. 

In a first step we present a detailed analysis of the electromagnetic form factors of 
the TT and a mesons. To solidify our model considerations, we study the electromagnetic 
form factor of the charged pion, which was recently measured by the Jefferson Lab 
Collaboration [43] . We also calculate the form factors which govern the transitions ^ 77 
with = 7r° and a including different kinematics of the photons. Our results for the 
pion are in good agreement with the recent experimental data [33,43]. We furthermore give 
results (analytical formulas and numerics) for the asymptotics of the tt ^ 77 and o" — > 77 
form factors at large values of space-like photon virtualities. The behavior of the 77^77* 
form factor nicely coincides with the prediction of perturbative QCD [32]. Expressions for 
the leading order (LO) ~ l/<5^ and the ncxt-to- leading (NLO) order ~ 1/Q^ expansion 
coefficients in the form factors are calculated at arbitrary values of large photon virtualities. 
We complete the analysis by indicating results for the charge radii. 

Based on the chiral symmetry construction we determine the a meson mass — 385.4 
MeV and the total width = 173 MeV. Both values are lower than the experimental 
values of the E791 [12] Collaboration, but close to the BES [13] results. The predicted 
strong coupling constant gf^TrTr = 1.8 GeV is fairly stable with respect to variations in the a 
mass, but the total width, due to phase space, depends sensitively on this value. Turning to 
the experimental evidence for the a meson we give predictions for the decay characteristics in 
the nonleptonic D ^ an and B an transitions. Both form factors and decays widths are 
determined. Our result for r{D~^ avr"'") = 0.90 x 10~^^ MeV is in agreement with the lower 
value of data [1,12]. Our prediction for the decay width T{B+ a7r+) is 0.94 x 10"^^ MeV, 
which is expected to be measured in forthcoming experiments. 
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APPENDIX 



A. Methods of calculation for matrix elements 



We demonstrate the evaluation of matrix elements for the example of the derivative of 
the meson mass operator (4): 



2p2 dp" J AttH 

The technique we use is based on the three main ingredients: 
• use of the Laplace transform of the vertex function 



(95) 



$(z) = Jds^Lis) e-'\ 





the CK-transform of the denominator 

1 



— {k + pY 



C50 



and the differential representation of the numerator 



After straightforward algebra wc get expressions for the coupling constants hn = l/n^(m^) 
of pseudoscalar {H = P) and scalar {H = S) mesons, which are given by 



oo 1 

, _i 1 f dtt f , r ,~2, T^f 1 ^1^2 , mj^ 



mim2 m\ A{1 — X^) + B 







(96) 



+ $^(.o)(l-A^ + ^)}, 

where [$^(2:0)]' = d[^'^{zQ)]/ dzQ, dn = for pseudoscalar mesons and dn 
mesons. We also define 



— 1 for scalar 



B. 



zo = j[(mi + m2f - m%] + ^(2a - 1 + A)^ 

A^l + t, B ^t[4:a{l- a)A + {2a-l + Af], A = ^21 - ^12 . 
In analogy we consider the basic vertex function (15): 
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Here we indicate the resulting expression for the Gaussian meson-quark vertex function 
t«(4)=exp(-4/A|,): 

;..3.3(p„p^) ^ - g„.) -P(-4)(|^ - . (98) 

where 

Ci = itr[(m2+ ^2)ri2(mi+ ^i)ri3(m3+ Pi)T2^] , 

C2 = 2^T['jTi2j„Ti3{m3+ p3)T23 + 7'^ri2(mi+ ^i)ri37^r23 + {m2+ -^2)ri27'^ri37^r23], 

Di^Pi-pt, D2^P2-Pt, Ds^-pt, Pt ^ ripi + r2P2 , 
ri = -^(toi + W13S13 + W12S12), r2 = -^{ta2 + W23S23 + W12S21) , 

3 

z = t[Y^ aimi - aipj - a2pl] + A[{ri + r2){ripj + r2pl) - rirsPg] 

i=l 

- Pli^ls^lS + W12S12) - ^2(^23^23 + W21S12) + P3W12W21S12 , 

1111 Ked 

— -I 1 s- — ^^'^ 



A2 A2 ■ A2 ' A2 ' •'^^ A2 ■ 

-'^-red ^^Hi2 ^^Hi3 ^^His ^Hj 

All further calculations are done by using computer programs written in FORM for the 
manipulations of Dirac matrices and in FORTRAN for the final numerical evaluations. 



B. Feynman rules for nonlocal electromagnetic vertices 

In the following we derive the Feynman rules for the nonlocal vertices of Figs. 3a and 3b. 
These vertices contain the path integral over the gauge field I{x, y, P). The crucial point is 
to calculate the expression 

^{dl)e^^[I{x,y,P)f (99) 

using Eq.(12). The case A'" = 1 corresponds to the vertex of Fig.3a and A" = 2 to the vertex 
of Fig. 3b. 

First, we consider the case A^ = 1. It is readily seen that 
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^dl)e'P^I(x, y, P) = e'^^^Vl)I{x, y, P) (100) 
where 'Dx = dx + ip- Thus we have to calculate 

Hvl)i{x, y,p) = Y: — y '^T n^, y, P) ■ (101) 

One finds that 

Vll{x, y, P) = [d:,A{x) + 2ipA{x)] - p^I{x, y, P) = L{A) - p^I{x, y, P) (102) 

where L{A) = d^A{x) + 2ipA{x). 
Iteration of the last expression 

{Viyi{x, y, P) = {Vl - p')L{A) + {-p'fl{x, y, P) (103) 
{Vlfl{x, y, P) = {Vl - Vy + p%{A) + {-p'fl{x, y, P) 

n-l 



{viri{x,y,p) = j: ivir-'-' i-p')' HA) + (V)"/(^,y,p) 

k=0 

1 

= n I dt[Vlt - p'{l-t)r-'L{A) + {-p'ri{x,y,P) 



finally leads to 

1 

^Vl)I{x,y,P)^ J dt^'[Vlt-p\l-t)]L{A) + ^-p')I{x,y,P) (104) 
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= / (^ ^'^^^^ {'^^^ + ^^""''^ / ^^^'["^^ + ^^'^ " ^'^^ " 

a; 

y 

where A^^q) is the Fourier-transform of the electromagnetic field and $'(2;) = d^{z)/dz. 
The last term in Eq.(104) containing an integration from y to x (path integral) vanishes due 
to the delta function S{x — y) in the Lagrangian. 

In complete analogy, we also obtain for the nonlocal vertex for the case N — 2: 

Hvi)i\x, y^n-J^J^, AM Mqi) (105) 

1 

X 1^29^" e^(«^+''2)^ J dt^'[-{p + qi + qift - p'(l - t)] 


1 1 

- (2p + qiY{2p + 2?i + q2y e^(^i+^2)- j dtt Ji 



1 1 

dl 
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X n-mp + p\i - 1)) - (1 - imp + q,f + (1 - t)p')] 

1 1 





X ^"[-l{t{p + qi + q2f + - t)) - (1 - + q^f + (1 - t)/)] 

1 a; 



+ i{2p + g2)'^e*^2" / citl>'[-(P + (12ft - /(I - t)] / dz'^e''^^ 





a; a; 

+ g^'^i-p") j dz^e''^'' J dw'e 
y y 



iq2W 



where ^"{z) = d'^^{z)/dz'^ . Again, the last three terms in Eq.(105) containing integrations 
from y to a; vanish due to the delta function. 



C. Form factors characterizing the local cr 77 diagram 

In the following we give the full analytical expressions for the form factors of Eq. (43) : 

(106) 



= ^2 {-2 - - ?2) A - 2Lig? (A - 6gi {ql - qf - p')) 

~2L,qU\-6ql {ql - ql - p")) 

+ 2 L,{ql + - q{ ql _ ql q^ - 2 q^p^ - 2 qtp^ + qlp^ + qlp^ + 8 ql qlp') 
+ hip', gl, ql) -ql- ql) (4m^ A - / + + ql - q\ ql - ql qi 
-3 gf - 3 q\ p" + 3 g2 p4 + 3 ql / + iq ^2 ^2 | ^ 



G£ = ^ {4 A + 4Li(A - 3g? (g^ + g^^ +/)) + 4L2 (A - 3g2^ {ql + g^^ +/)) 

+4L3(A-3/ {ql + ql+p')) 
- 2 I^{p\ ql ql) (4 A + / + g^ + ql - qt ql - ql qt -qtp'- qiv" 

-qlp'^ -qlp'^ + Qqlqlp'^)} 

where /l(p^, ql, ql) is the integral defined in Eq. (37). Here we use 

1 



Li = j dx In 



with g3 = p. 



1 — x{l — x) 



= -2 + 2 



4 — g| 



ArcTan 



4 777,2 _ q2 



(107) 



(108) 



D. Two-loop integral representation and power expansion of /nl 

The integral /nl of Eq. (46), which is related to the 77^77 form factor, can be reduced 
to the two-loop integral 
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J l + tJ \tAi(x) + Ai tA2{x) + Ai tA3(x) + Aif' ^ ' 



where 



Dp — rri^ — X {1 — x) 
Dg. =m'^-x{l-x) qf 
Ao = + pV4 - {ql + ql)/2 
Ai = -pV2 + (g? + ql)l2 
A2{x) = (1 - ao) 2; • - (1 - Qio) gi - cto gi 
^3(2^) = -ao a; • A^/^ - (1 - ao) gf - olq ql 

With the constraint that < 4m^ and qf < the variable t is replaced by m = 
+ t) +t A/(l + t) in the integral of Eq. (109): 

00 1 \/-\7-)/-\ 

Jnl = Idxi \- . ^" ~ 1/5; ! w- (110) 

I /I i?j,(l) [ii-A(l) + i?p(l)]-Ai(x) + 2DpAi ^ ^ 

(1-Qi0)[^- A(x)+i^g,(x)] 

i?,,(a:) [ti - A{x) + i?g,(x)] • A2{x) + 2 L>g, Ai 

1 q;o[-»- A(x) + i?g,(a;)] ) 

i?,, (x) [li - A(x) + Rq, {x)] ■ A^ix) + 2Dg,Aj' 

where Rp{x) = {u - A{x)Y + AuDp. 

Next we perform a power expansion of /nl in the limit — For this kinematics we 
have 

D = + \az Q2 (ai (1 + a;) + q;2 (1 - u)) , (111) 
A = + • 
By scaling the variable 0:3 — > a^/Q^ one finds 

- -^/* (1^)7"'"^ - - - - - i) *' (""^ + ^) • ("^' 

where 

W ^l + t[ai{l + uj) + a2{l-uj)]. (113) 

Prom the last expression it is easy to determine the asymptotics in the leading order (LO) 
with 



zp^ ^,Dp+^^ A{1) 

A(a;) = Ao + a;Ai 
Ai{x) = x ■ A^/^ + ql - CioP^ 
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(^nl)lo = (114) 



By using the formulas 

j(fa d'{l — ai — 0:2) F{ai, 02) = j(fa 6(1 — ai — 0:2) (-^(0, 1) + F^^ {ai, a2 

(115) 

00 CXD 00 t 



dtg{t) j dr f{T) = Jdtf{tm^) jdrgir), 



tw? 00 

where F^^(q;i,q;2) = dF{ai,a2)/dai, the next-to-leading (NLO) term is obtained as 

Note that both integrals B}f{oj) and i?J^^°(a;) are defined in Eq. (51). 

E. Form factors characterizing the gauge-invariant part of the nonlocal triangle 

cr 77 diagram 



Here we give the full analytical expressions for the form factors of Eqs. (61) and (62): 

m 
A2 



^Xx = 5 - + ^ + 2 ql (A - 12 + 12 g^'p^) (117) 



-2Lf (g? + ql - qUl " -2qtp' -2qtp' + g?J>^ + qlp' + SqlqW) 
+6 1^11 {ql + ^2' - P') {<il + ^4- P') - 6 V^2 ql {ql - ql - p'M + ql " p') 
+6 1^21 ql {ql + ql - p') (3 ql + ql - /) + 6 V22 ql {ql + ql - p'M - ql + p^) 
+/NL(P^ ql ql) (/ - ql - ql) (4 A - / + + - gf g,^ - g^ g^ 

-3 g^ p2 _ 3 ^4 ^2 ^ 3 ^2 p4 ^ 3 ^2 ^4 ^ ^2 ^2 ^2 ^ | 



^A, - 1-^ {S, - S,) {ql + g2^ - /)(A + 12 ql ql) (118) 
91 92 



m 

Pi; 

+4 L^'^(g^ - 2 g,^ - 2/ - 5 ql ql + g?/ + 4 g,^/) 
-4 (2 g^ - g,^ + 2 + 5 ql ql - 4 g,^ - gi p^) 
-4 (g^ + g,^ - 2 / - 2 g? g,^ + g? / + ql /) 

+2 Vn^(g2 + 3 g2 - pl{\ + 6 g^ g^) - 2 V^i2^(gi' - g2' - /)(A + 6 g^ ql) 

92 92 

+2 F2i^(3 ql + gi - p^)(A + 6 g? g^) + 2 1/22^(9? - 92' + /)(A + 6 g? g^) 

-2 J^L(k 9?, 9i) (4m^ A + / + g,^ + ql - qt ql - ql qt - qtp' - q^p' 



2 4 2 4,^- 2 2 2\1 

-9iP -92P +6gig2P )] 
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where /nl(p^, ^i) is the three-dimensional integral of Eq. (46). Above we introduce the 
terms 



rNL 



oo 1 



^ ^ 



oo 2 1 



= Jdt-^^—r^ jdxx ^{Zi) , (121) 



^ ^ 



oo 2 1 



= j y (1 - ^) ^(^0 , (122) 

^ 



where 



t 

/ 2 /I \ 2\ , ^ ( 2 , t 2 



^(m^-x(l-x)g,^) + ^p-^ + |(p^-g?-g2^)j , (i = l,2), 



Z3^Y^{m -x{l-x)p ) + j-^\rn - — 



The asymptotics of these expressions can be derived from their representation in a manner 
as discussed for the 7r77 form factor in Appendix D. Here we obtain 



,3 



K).co=-il«-°(-). 

The integrals R^l'^{uj) and i?^^*^(co') are defined in Eq. (69). The last simple relation between 
(^Faj^) lo ^'^'^ ('^^J^) Lo <i6^ived with the help of the identity 

tu{ai - a2) = W - {1 + t («i + 0^2)), (124) 
where W is defined in Eq. (113), and the exchange symmetry cti q;2. 



F. Form factors characterizing the gauge- invariant part of the nonlocal bubble and 

tadpole (T77 diagrams 

The gauge invariant parts of the nonlocal (777 bubble (bub) diagram introduced in Eq. 
(59) are written as: 
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OO r, 1 1 

m r , r 



= y / ^^71-7)4 Jdrrjdxil -2x){^'{z^) + $'(^2)) , (125) 



[1 + t) 

^ ^0 

00 11 
m r t 



1 



/ /"^^ /"^^ {^^'"'^ ^ ^^'"'^ ^^^^^ 







+ 4 



iv'-tl- 92)^^(1 - 2i)(4'fe) + i'fe))} , 



where 



_ T(l-r + t) g| _ ^ 2 9? 

1 + t 4 1 + t ^ ^ 4 ' 

Zi = z{qi,q2), Z2 = z{q2,qi) . 

The final expressions for the tadpole (tad) diagram, originally defined in Eq. (60) , are 

00 11 

^ ^0 

X 



||.'(Z0 + ^'{Z2) - ^'{Zs) - ^'{Z,)] (1 - -^^j , 

00 1 1 



^ '0 

^2 



+ |i+^^)2 ^^ (gi + ^2 -^^') + *'(^2) - $'(^3) - $'(^4)) 

^ ($(zi)+l'(z2)-$(z3)+$'(^4)) 



1+t 

where 

2 p^l + t — r 2 TTr7-.l + ^~''' 

^(p) = m t - — T , Zp = z{p) , 2;^ = z{q) , 2;^ = m t - PFi - i^i , 

W(g) = ^ r / (1 - /) , i?(gi, q2, 1) = ^ {l [p' - ql - ql] + ql + qt) , 

Ri = R{qi,q2,l) , R2 = Riq2, qi, I) , R3 = R{qi, q2, -I) , R^ = R{q2, qi, -I) ■ 
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TABLES 

TABLE I. Leptonic decay constants fp (MeV) used in the least-square fit. Data are either 
taken from PDG [1] or from the Lattice [27] (quenched (upper fine) and unquenched (lower line)). 



Meson 


This model 


Expt/Lattice 


7r+ 


131 


130.7 ±0.1 ±0.36 


K+ 


161 


159.8 ±1.4 ±0.44 


D+ 


211 


203± 14 






226± 15 


Dt 


222 


230± 14 






250± 30 


B+ 


180 


173± 23 






198± 30 




196 


200± 20 






230± 30 


Bt 


398 
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FIGURES 



Fig.l: Meson mass operator. 

Fig. 2: Hadronic decay H12 H12 + H23. 

Fig. 3: Nonlocal coupling of meson, quark and photon fields: 
one-photon vertex (3a), two- photon vertex (3b). 

Fig.4: Diagrams contributing to the electromagnetic form factor Ft^{Q'^) of the pion: 
triangle (4a) and bubble diagrams (4b and 4c). 

Fig. 5: Diagrams contributing to the processes H 7*7*: 
triangle (5a), bubble (5b and 5c) and tadpole diagrams (5d). 

Fig. 6: Diagrams contributing to the decay amplitude D[B)^ utt"'": 
direct (6a) and polar (resonance) contributions (6b). 

Fig. 7: Electromagnetic form factor Ft^{Q'^) of the pion in comparison to data taken from [43] 
(JLAB), [46] (DESY), [47] (CERN) and to other theoretical calculations: QCD sum rules 
[48], light-cone quark model [49], NJL model with a separable qq interaction [50] and QCD 
modeling approach based on the solution of the Dyson-Schwinger equations [51]. 

Fig. 8: Results for the form factor Q"^ F-,^^^* {Q"^) in comparison to experimental data taken 
from [53] (CELLO), [33] (CLEO) and to other theoretical calculations: hard scattering 
approach (HSA) [54], QCD sum rules [55] and perturbative hght-cone QCD [56]. 

Fig. 9: Our result for the form factor Q'^F^^*^*{Q^,u!) for different values of the asymmetry 
parameter uj. 

Fig. 10: Our result for the form factor Q"^ F^^jj* {Q"^) . We also indicate the separate contribu- 
tion of the dominant triangle (A) diagram. 

Fig. 11: Our result for the form factor Q'^F^^*^*{Q'^ ,uj) for different values of the asymmetry 
parameter uj. 
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